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Abstract 

We apply the Linear Delta Expansion (LDE) to the Lindstedt-Poincare ("distorted time") 
method to find improved approximate solutions to nonlinear problems. We find that our method 
works very well for a wide range of parameters in the case of the anharmonic oscillator (Duffing 
equation) and of the non-linear pendulum. The approximate solutions found with this method are 
better behaved and converge more rapidly to the exact ones than in the simple Lindstedt-Poincare 
method. 
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I. INTRODUCTION 



The study of nonlinear problems is of crucial importance in all areas of Physics. Some 
of the most interesting features of physical systems are hidden in their nonlinear behavior, 
and can only be studied with appropriate methods designed to tackle nonlinear problems. 
In general, given the nature of nonlinear phenomena, the approximation methods can only 
be applied within certain ranges of the physical parameters and or to only certain classes 
of problems. It is a challenge to devise nonlinear frameworks that contain both operational 
ease and flexibility in their application. In this paper we present a method for the solution 
of nonlinear problems that attempts to accomplish these features. 

There are several methods which have been used to find approximate solutions to nonlin- 
ear problems. Here we just review a few. Lindstedt developed a method long time ago Q in 
which one considers solutions to problems involving conservative oscillatory systems with an 
unknown period. The main observation is that by introducing a rescaled time, one can avoid 
the appearance of terms indefinitely growing with time ("secular terms"), that are common 
in ordinary perturbation theory. The method is now know as the Lindstedt-Poincare (LP) 
method or as the Distorted Time method. 

Another known technique is the perturbative 5 expansion (see for example Q). In this 
case the idea is to modify the exponent of the nonlinear term by introducing a parameter 
5 as new exponent. 5 interpolates between the linear (5 = 0) and the nonlinear (5 = 1) 
problems. If one is able to solve the linear problem then the original nonlinear problem 
becomes, after a power expansion in 5, an infinite sequence of linear problems which are 
(formally) solvable. 

n 

Yet another framework is the Multiple-Scale Perturbation Theory (MSPT) [3]. In this 

case, one tackles problems in which a dynamical system has physical behaviors at various 

length or time scales. This is usually problematic for ordinary perturbation theory due to 

the appearance (again) of secular terms. The central idea is to introduce more than one 

time and to treat them as independent variables. By performing the usual perturbative 

expansion, one then imposes conditions on the solutions (which depend on the different 

"times" ) in order to get rid of secular terms and a linear differential equation is left to solve. 

fl 

Finally, the Linear Delta Expansion (LDE) |4j|. This is a method in which an arbitrary 
(or several) parameter A is introduced into the problem and calculations are carried out 



with conventional perturbation theory in an expansion parameter 5=1. At each order in 5, 
the convergence of the approximation can be improved by applying the principal of minimal 
sensitivity which consists on a minimization of an observable with respect to the parameter 
A. 

n 

All of these methods have been applied to a variety of problems. In |2[, Bender et al. 
showed how one can obtain approximate solutions using the perturbative S expansion and 
the MSPT to the Duffing equation (the classical anharmonic oscillator). Its success then 
has motivated their extension of the method into quantum systems jjjj. The LDE method 
has extensively been applied in many different settings with varying degrees of success. For 



rely 

Eh 



example, in 51 it has been used to analyze disordered systems. Pinto and collaborators have 

□ 

applied it to the Bose-Einstein condensation problem p, the 0(N)(<f) 2 )1 d model jjjj, to the 
Walecka model Q and to the <f) A theory at high temperature Js|. Detailed references can be 
found in these works. 

We can see that it is possible to tackle a large number of nonlinear problems with these 
well known techniques. However, there is still room for substantial improvement over them. 
As mentioned before, it is desirable to have a method that works over a large range of 
parameters, which is not always the case in the aforementioned methods, and we would like 
the new method to give a smaller error in the approximations than its competitors. It is also 
desirable to devise a framework with operational flexibility and so easy to adapt to many 
different problems. 

We show that the method presented in this paper accomplishes these features in the 
case of the Duffing equation and of the nonlinear pendulum. The method is based on the 
application of the LDE to the LP method jl^]. We find solutions that are better behaved 
and that converge much faster than in the other methods described. 

In Section |H] a brief review of the LP method is presented followed by a review of the 
LDE method in Section ITTT1 We then show the application of both methods to two problems, 
the Anharmonic Oscillator in Section IIVI and the nonlinear pendulum in Section We 
present our conclusions and current work in Section IVII Appendix A contains some of the 
formulae employed in the computations. 



II. THE LINDSTEDT-POINCARE METHOD 



a. 



In this section we introduce the Lindstedt-Poincare distorted time (LP) method [1]. We 
consider a nonlinear ODE of the form 



x(t) +uj 2 x(t) = e f(x(t)) 



(1) 



which describes a conservative system, oscillating with an unknown period T. The nonlinear 
term e f(x(t)) is treated as a perturbation. Unfortunately, when the ordinary perturbation 
is applied to eq. ((TJ), by writing the solution as a series in e, the appearance of secular terms 
spoils the expansion and any predictive power is lost for sufficiently large time scales. 

In order to avoid the appearance of secular terms, we switch to a scaled time r = 2irt/T = 
Q t, where T is the (unknown) period of the oscillations. The ODE now reads: 



fi 2 ^(r)+a; 2 x(r)=ef(x(r)). 



(2) 



We notice that the dependence upon e in this equation enters both in the solution x{r) 
and in the frequency Q. By assuming e to be a small parameter we write 



tt 2 = £ e n a n ; x{r) = £ e n x n (r) 

n=0 n=0 

and expand the r.h.s of eq. (J2J as 

/(*) = / ( E z n x n(r)) « f(x ) + e Xl f(x ) + e 2 



\n=0 



X2 f'(x ) + -f f"(x ) 



+ s 6 



X3 f{x ) + X 2 Xi f"{x ) + -± f'"(x ) 





+ 



By using these expansions inside eq. (j2J) we obtain a system of linear inhomogeneous 
differential equations, each corresponding to a different order in e. Let us consider the first 
few terms. To order e° we obtain the equation 

d 2 x 



Q'o 



dr 2 



describing a harmonic oscillator of frequency fl 
equation 

d 2 x\ 



+ oo 2 x (t) = , (3) 
= to. To order e we obtain the 



a 



dr 2 



+ u xAt) = Si{T) 



(4) 



where the r.h.s. is given by 

si(r) = -ax -^y + f(x ) . (5) 

We stress the oscillatory behavior of the driving term Si(r), because of its dependence 
upon the order-0 solution, Xq(t). As a result si(r) will contain the fundamental frequency, 
corresponding to a period of 2ir in the scaled time, and multiples of this frequency, appearing 
through the term f(xo(r)). The presence of a driving term with the fundamental frequency 
leads to a resonant behavior of Xi(t) and to the unfortunate occurrence of secular terms, 
which spoils our expansion. However, we can deal with this problem by fixing the coefficient 
oi\ to cancel the resonant term in the r.h.s. of eq. The iteration of this procedure to 
a given order n allows to determine the coefficients ao, . . . ,a n and therefore the frequency 
Q = y/a + e «! + ... + e n a n . 



III. LINEAR DELTA EXPANSION 



The linear delta expansion (LDE) is a powerful technique which has been originally 
introduced to deal with problems of strong coupling Quantum Field Theory, for which the 
naive perturbative approach is not useful. Since then this method has been applied to a 
wide class of problems 0, Q, 0, 0]. In its original formulation a lagrangian density C, 
which is not exactly solvable, is interpolated with a solvable lagrangian £o(A), depending 
upon one (or more) parameters A: 

Cs = C (\)+5 (£-A,(A)) . (6) 

For 5 = one obtains Cq(\), whereas for 5 = 1 one recovers the full lagrangian £$. The 
term 5 (£ — Co) is treated as a perturbation and 5 is used to keep track of the perturbative 
order. Eventually 5 is set to be 1. 

We notice that the interpolation of the full lagrangian with the solvable one, £o(\), brings 

an artificial dependence upon the arbitrary parameter A. Such dependence, which would 

vanish if all perturbative orders were calculated, can be soften to a finite perturbative order, 

by requiring some physical observable O to be locally insensitive to A, i.e: 

dO(X) _ 
d\ 

This condition is known as Principle of Minimal Sensitivity (PMS) and is normally seen to 
improve the convergence to the exact solution. 



IV. ANHARMONIC OSCILLATOR 



In this Section we apply the LDE to the LP method in order to find approximate solutions 
to the Duffing equation, a problem which has already been considered in [lOj]; here we present 
the calculation in more detail. 

Consider the equation for the anharmonic oscillator 



-fi x 3 (t) . 



This equation describes a conservative system, where the total energy is given by 



E= X - + 
2 



to 2 x 2 



X 

2 + /i T 



The period of the oscillation can be calculated in terms of an elliptic integral 



exact 



dx 



y/2(E-V(x)) ' 

where A is the amplitude of the oscillations. 

Following the procedure explained in the Section |H] and IHIl we write Eq. (J7J) as 
,d 2 x 



(?) 



(8) 



(9) 



dr 2 



(r) + (oo 2 + A 2 ) x(r) = 5 [-/j x 3 (t) + A 2 x(r) 



(10) 



where an arbitrary parameter A with dimension of frequency has been introduced. Clearly 
for 6 = 1, Eq. (fTUj) reduces to Eq. ((7j). We repeat the procedures previously explained and 
find a hierarchy of linear inhomogeneous differential equations to be solved sequentially. 

1. Zeroth Order 



To zeroth order we obtain the equation 



with solution 



Xq(t) = A COST 



(11) 



(12) 



The zeroth order frequency is then given by 



a = oj 2 + A 2 



(13) 



2. First Order 



To first order we find the equation 



where 



S'i(r) = A cosr 



«i + A 2 



3 Ay 



cos3r . 



Now ai is fixed by eliminating the term proportional to cos r: 



A 2 . 



We obtain the solution 



. . A 3 n A 3 jj, 

Xl(T) = - 32(o;' + A*) COST+ 32RTP) COS3r ' 



and the frequency 



Q 2 = a + ai = u 2 + 



3A 2 fi 



which is observed to be independent of A. 



3. Second Order 



The second order equation is given by 



where now 



«o ^ + (^ 2 + A 2 Wt) = S 2 (t), 



. . A (3 A 4 jj, 2 + 128 a 2 (uj 2 + A 2 )) 

5a(r) = 128 p + A') COST 

A 3 /x(3 A 2 /i-4A 2 ) 

+ i fi / 2 i \2^> C ° S3r 

16 (or + A 2 ) 
3 A 5 fi 2 

cos or . 



128 (lu 2 + A 2 ) 

As before a 2 is fixed by eliminating the term proportional to cosr: 

/i 2 



3 A A •■ 2 



a 2 = 



128 {oo 2 + A 2 ) ' 



We obtain the solution 



A 3 fi (23A 2 n - 32A 2 ) A 3 fi(-3A 2 fi + 4A 2 ) 

X2(T~) = ttttt: — COS T H -— — — - — COS OT 

A ; 1024(cj 2 + A 2 ) 2 128(c; 2 + A 2 ) 2 



+ Wm$TW- cos5t (21) 



and the frequency 



, , 3A 2 fi 3A 4 ^ 2 



n=a » + ai + a2 = w+ — -128FTP)' (22) 

Note that at this order the frequency now depends on the arbitrary parameter A. However, 
due to the explicit dependence, by applying the PMS, we would obtain the same solution 
as in the simple LP method (A = 0). In order to get a different solution, we must go to the 
next order in the expansion. 



4- Third Order 



Following the same procedure, we obtain the following expression for the third order: 



"0 + (" 2 + A 2 ) s s (r) = S 3 (r) 



where 



ss(r) 



COST 



3 A 5 /1 2 (3 A 2 n - 4 A 2 

A a% — 

L 3 512 (w 2 + A 2 ) 2 

(A 3 fi (297 A 4 /1 2 - 768 A 2 /1 A 2 + 512 A 4 )) 

7—z — cos3r 

2048 {uo 2 + A 2 ) 2 



3 A 5 fi 2 (3 A 2 fj, - 4 A 2 ) 

+ 7—7; — cos5r 

256 (A 2 + w 2 ) 2 



3 A 7 /i 3 



POS ) — — . . ... - 

2048 (A 2 + uo 2 ) 2 
we determine 0:3 to be 



cos 7r . 



By eliminating the term proportional to cos r 

3 A 4 /i 2 (3 A 2 /I - 4 A 2 ) 

«3 — 



512 {\ 2 + uo 2 f 
and the solution 

A 3 fi 547 A 4 fi 2 - 1472 A 2 n A 2 + 1024 A 4 
Xs(t) = "32768 l^TU 2 ) 3 COST 



+ 



A 3 n 297 A 4 fi 2 - 768 A 2 [i \ 2 + 512 A 4 
T7^7 (A 2 + ^ 2 )3 COS3r 



16384 v „ ^ , 

A 5 fi 2 (-3 A 2 fj, + 4 A 2 ) A 7 /i 3 

+ 2048 (A 2 + uo 2 ) 3 C ° S + 32768 (A 2 + co 2 f 



cos 7r . 



(23) 



(24) 



(25) 
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FIG. 1: Squared frequency of the anharmonic oscillator as a function of the amplitude (arbitrary 
units), u = fj, = 1. 



The frequency to order S 3 is now obtained to be 

3A 2 /j 3 A A fi 2 



| 3iV 2 (3AV-4A^ 



4 128 (cj 2 + A 2 ) 512 (A 2 + cj 2 ) 2 

This time, the frequency depends upon the arbitrary parameter A in a nontrivial way and 



we can apply the PMS in order to fix the value of A. We do this by imposing that 
which leads to the following result: 

„ 3 A 2 a 
A = — — - . 



dX 



0. 



(27) 



Notice that since A depends linearly upon A the formula for Q 2 obtained in this case does 
not simply correspond to an expansion in A. As a matter of fact we find that the frequency 
corresponding to this value of A is 

64 A 4 /i 2 + 192 A 2 fi uo 2 + 128 u A 



(28) 



96 A 2 fi + 128 uo 2 

Notice that the Duffing equation (0) is left invariant under the simultaneous rescaling of 



the anharmonic coupling fi and of the amplitude, i.e. and A — > A' = A y ft/ ft'. This 

invariance is manifest in the equation (|2*H|) . which is function of A 2 ft, which is invariant 
under this rescaling. 
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FIG. 2: Squared frequency of the anharmonic oscillator as a function of the amplitude (arbitrary 
units), lu = 1 and fj, = — 1. 

In Fig^ we compare the exact frequency, calculated with Eq. (JHJ) with the frequency 
obtained with our method (LPLDE), equation (|28jl. and with the LP method, equation (J2HJ) 
taking A = 0, both to third order in perturbation theory. We take u = \x = 1 (see the left 
plot of Fig. EJ) and vary the amplitude of the oscillations. We observe that our method yields 
an excellent approximation to the exact result even for large amplitudes, where the simple 
LP approximation fails. 

In FigElwe consider the case studied in FigJU but choosing u — 1 and /i = — 1 (see the 
right plot of Fig. In this case the potential has a local minimum in the origin and two 
maxima, located at x — ±1. Periodic solutions are supported only for amplitudes A < 1, 
A = 1 being a point of (unstable) equilibrium, where the period diverges. Also in this case, 
the LPLDE method offers an excellent approximation to the exact result for a large range of 
amplitudes; as expected, the approximation is poorer in the region Aral, where the point 
of equilibrium is approached. 

In Fig0]we compare the period obtained with our method to the exact period of Eq. (jHJ) 
and to the one obtained with the formulae of Q], which are obtained by applying the nonlin- 
ear delta expansion. Our method provides an excellent approximation to the exact period 
over a wide range of the parameter /i, which controls the nonlinearity. The plots are ob- 
tained assuming u = 1 and the boundary conditions x(0) = 1 and x(0) = 0. The formulae 
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FIG. 3: Anharmonic potential corresponding to A) u = fj, = 1 and B) u = 1 and = —1. 

of j^J behave badly in the region /i < 0, which corresponds to a potential well of finite 
depth centered around x = 0, and yield a precision comparable to the one achieved with 
our method for fi > 0. Corresponding to the value \i = the oscillator is in a position of 
(unstable) equilibrium and the exact period diverges. Notice that for large values of /i all 
the methods seem to give a good approximation to the exact solution, including the LP 
method (to first order), which (to third order) was behaving poorly in the case previously 
studied. Unfortunately the equations of Q are not suitable to be analyzed as in Fig. [TJ and 
thus we cannot fully test the efficiency of this method. 

In FigElwe plot the relative error corresponding to the different approximations for > 0. 
Our method to third order in perturbation theory yields an error typically smaller than the 
errors of the other methods and with a magnitude of about 0.1 %. 



V. THE NONLINEAR PENDULUM 



We now apply the improved method to the nonlinear pendulum. The steps are exactly 
the same as before and we proceed to outline them. First, consider the equation for the 
nonlinear pendulum 

d 2 9 

— + uj 2 sinO = 0, (29) 
at 2 

where uo 2 = g/l is the natural frequency of the small oscillations of the pendulum. Following 
the Lindstedt-Poincare method, we introduce a scaled time r = Q t and write the equation 




Exact 

BMPS (order 1) 
BMPS (order 2) 
LP (order 1) 
LPLDE (order 3) 



(H+l) 



FIG. 4: Period of the anharmonic oscillator. The curves labeled with BMPS refer to the formulas 
of 



Vjr. 

El 



as 



Q 2 — + uj 2 sin 9 = 
dr 2 



(30) 



where Q = 2ir/T is the (unknown) frequency and T is the period of the oscillations. As 
discussed in the case of the anharmonic oscillator, we can apply the Linear Delta Expansion 
to the problem by modifying the above equation and writing it as: 

d 2 e 



dr 



] +\ 2 e = 5 -u 2 sm6 + \ 2 6 =Sf{9) 



(31) 



where A is an arbitrary parameter, with the dimension of frequency. In what follows we use 
the same procedure previously outlined for the anharmonic oscillator, with a few technical 
differences due to the more difficult nature of the present problem. 
We will expand the angle and the frequency as 



e{r) = E r W 



n 2 



n=0 



n=0 



We will solve Eq. (|31[) subject to the boundary condition 9(0) = A and 9(0) = 0, i.e. 



9 (0) = A , 9 J>0 (0) = , 9 j = 



(32) 
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FIG. 5: Error corresponding to the different approaches for the case studied in Fig. 21 The curves 
labeled with BMPS refer to the formulas of 2|. 



5. Zeroth Order 

To zeroth order the equation for the pendulum reads 

«o + A 2 9 = , 

and we obtain the solution 



(33) 



Oq(t) = A COS T 



(34) 



describing a simple oscillatory motion with (scaled) period 2ir. The zeroth order frequency 
is therefore given by 



a = A 2 



(35) 



6. First Order 



To first order we obtain the differential equation: 



"o^ + A 2 ^ = S 1 (t), (36) 



where we have defined the source term 
d 2 6 



Si(r) 



-Oil 



dr 2 



f(6 ) = Aai cost+ —u sin (A cost) + A A cost . (37) 



As before, in order to avoid the occurrence of secular terms, we need to eliminate con- 
tributions proportional to cost from the source term Si(t) (recall that such a term would 
yield a resonant behavior of the solution 9i(t)). We enforce this condition by requiring that 



1 r 2n 

- / dr Si(t) e tT = . 

7T JO 



As a result of this operation, we are able to fix the coefficient ai~. 



ai = A 



00' 



2tt 



-A A H / dr sin(A cost) e l T 

7T JO 



1 

A 



-A 2 A + oo 2 a 



where we have used the following expansion of sin(v4 cost): 



and 



sm(A cost) = E c 2j+i cos [(2j + 1)t] 

3=0 



c 2j+ i = - dr (2i+1)r sin(A cost) = 2 J 2j +i{A) 

TV JO 



Eq. (|3Tj) now reads 

d 2 6i 
~d^ 



Oi 



+ \ 2 6i = Si(t) 



-oj 



Yl c 2i+i cos (( 2 i + 1 ) T ] 



(38) 



(39) 



(40) 



(41) 



(42) 



where the sum starts from j = 1 because of the vanishing of the term proportional to cosr. 
We write the solution #i(t) as 

oo 

di(r) = E 4 } + i cos[(2j + l)r] , (43) 

3=0 



where the coefficients are (for j > 1): 



(1) _ u c 2j+ i 
23+1 ~ 4 \*j (j + 1) 



d 



(i) 

2j+l 



A 2 



(44) 



In the last equation we have introduced the scale coefficients dfj +1 , which do not depend 
upon A. We notice that Eq. ()42)1 cannot be used to determine the coefficient corresponding 
to j = 0; in fact, this coefficient is fixed by the boundary condition: 

oo 

7(1) _n (45) 



0i (o) = E 4+i = o 

3=0 



which entails 



d\ 



(i) 



V d (l) 

2^ a 2j+l 



E 



OJ C 2 j + l 



d\ 



(1) 



i=i 



{ 4 A 2 j (j + 1) A 2 



(46) 



7. Second Order 



To second order we obtain the equation: 



d 2 9 



«o ^ + A 2 9 2 = S 2 (t) , (47) 



where we have introduced the source term 

S 2 (r) = -«i§-« 2 §^ l(T) /W 

= "«i - «2 ^ + 0i (r) [-co 2 cos 0o (r) + A 2 ] . (48) 
We can expand the source term in a series as 

oo 

S 2 (t) = E 42+1 cos(2n+l)r , (49) 

n=l 

where the coefficients of the expansion are given by 

1 2tt -( 2a ) 

42 + i = - / (-+^ 5 2 (r) = + 4Si • (50) 

7T JO A Z 

We have introduced the scaled coefficients s 2 n+i an d 4n+i> which are independent of A and 



read: 



«i = ^(2n + l) 2 4 ) + i 



— \ 2^ a 2j+l C 2(j-n) + 2^ «2(/-n-l)+l C ^ + 2^ "2j+l C 2(n-j) ( , 
j=n l=n+l j=0 



2 

s 2 nji = -4 n {n + 1) rfsn+l • 
The coefficients c 2 j follow from the expansion of cos [A cos r] : 

oo 

cos [A cos t] = E %• cos [2j r] (51) 

j=0 

and read, for j > 0, 

1 / >27r 

c 2 j = — dr cos [A cos r] e* 2j r 
7T Jo 

2E(-ir( ¥ ) J— J, =2 { -iy M A) (52) 



and, for j = 0, 



c = cos A — E c 2j • (53) 

3=1 



(2) 

As before we need to eliminate the coefficient s\ 



2 oo 2 

a x d? + a 2 A + \ 2 d? 4+1 (%' + C2 j+ 2) ~ %r d? 5 (54) 



2 1 1 

3=0 



and obtain the coefficient a 2 



« 2 = 34 [ - r)lF + 2X^ d ^ {C2 > + Cw) J " P ' (55) 

where «2 = A 2 a2 is a scaled coefficient, independent of A. 
We are therefore able to find the solution of Eq. (|47p 

oo 

#2(t) = £ 4?+i cos [(2j + l)r] , (56) 

3=0 



with the coefficients, for j ^ 1 



j(2«) j( 2b ) 
j(2) _ »2j+l a 2j+l 

expressed in terms of the A-independent terms: 

-(2a) 

-pa) _ _ S 2j+1 

2J+1 ~ 4 j (j + 1) ' 

-(26) 

-,(26) _ g 2j+l _ t(1) 

S+l " 4 ^ (j + 1} " S+l ' 

As before the j = coefficient is not fixed by the equation and needs to be determined 
by enforcing the boundary condition #2(0) = 0. We obtain: 

00 00 „(2) 

^- r ? C 1w <58) 



8. Third Order 



To third order we obtain the equation 

d 2 9 3 



aiQ 



dr 2 



+ A 2 9 3 



(59) 



where the source term 63 (t) is 
S 3 (t) = -an 



d 2 9> 



dr 2 



d l B x d 2 9 , 
— «2 —rr, «3 —rrr + 



dr 2 



dr 2 



2 (r) f\9 ) 



f(0 ) 



(60) 



Once again it is useful to expand the source term in a series as 

oo 

Ss(r) = £ 41i cos(2n+l)r, (61) 

n=0 

where the coefficients of the expansion are given by 

, 2qT _(3o) -(3b) 

41i = ~ I dr j ^ S 3 (r) = ^ + ^ + 4Si (62) 

and s( 3a,fe,c ) are independent of A. A lengthy calculation allows to find the expressions for 
these coefficients, which can be found in Appendix A. Here we only write the coefficient of 
the term cost, corresponding to n = 0: 

up' 00 cu 2 00 

s { 3) = ai 4 2) + "2 4 1} + «3 -4 - — £ c 2j 4?+i ;r~ zZ ^2/ 4?-i 

2 j=0 2 1=1 

- ^ co 4 2) + a 2 4 2) 

, ,2 oo oo ,2 oo 3 

"7T Z^ c 2(j-m-l)+l «2m+l «2j+l + "7T Z^ Z^ C 2(-m+j)+l «2m+l «2j+l 

8 m=0j=m+l 8 j=0m=0 

, ,2 oo oo , ,2 oo m 

— 2^ c 2(m-j-l)+l «2j+l U 2m+1 ~>~ ~TT 2-^2-^ C 2(m-j)+l «2j+l a 2m+l 

8 m=j+l j=0 8 m=0 j=0 

^Jl oo oo ^J2 oo oo 

+ "7T X ZZ C 2(m+i)+l ^2j+l ^2m+l + "7T X X] C 2(m.+j+l)+l ^2j+l ^2m+l ■ (63) 
8 m=0 j=0 8 m=0 j=0 

The coefficient a% is fixed by requiring that sf^ vanish: 



+ 
+ 



"3 = + ^ ' (6 4 ) 



where 

"3 



UJ 2 j/Ci C \ -j(2o) _ -j(l) If,. ,. x -j(2a) 
«a = \\~A~~2J 1 + "2 ^1 - ^ Z^ ( C 2i + C 2i+2) rf 2j+ i 

1 °° °° (1) (1) 

+ o ZZ ZZ ( 2 C 2(i-m-l)+l + C 2 ( m +j)+l + C 2 (m+j)+z) d 2j+1 C^m+l 

8 m =0 \j= m +l 

m 

+ E ( 2 C2(m-j)+l + C2(m+j)+l + C2(m+j)+z) G? 2 j+1 ^2m+l 
3=0 

\ /Ci C \ -,(26) If -(26) \ _ , . 

aa = _ i IriJ 1 ~ 2 ^ ^ C2i C2j+2 ' ) 2j+i j = a2 ■ ^ ^ 

To this order the squared frequency reads: 

Q 2 „ , o «2 , «3a 

ft = aia + 2 — + — . 
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FIG. 6: Period of the nonlinear pendulum as a function of the amplitude. We assume to = \fgjl = 1. 
The "principle of minimal sensitivity" yields the solution 

A 2 = (66) 

and a corresponding value of Q 2 : 

= a la - . (67) 

In Fig. IHlwe plot the period of the nonlinear pendulum as a function of the amplitude, 
as obtained in the LPLDE and LP approximations, and compare the results with the exact 
period. We assume u — 1 and use the formulae given above truncating the infinite series to a 
maximum value j max = 5. As it can be seen from the Figure, the LPLDE approximation is in 
excellent agreement with the exact result, up to very large amplitudes. A = ±7r corresponds 
to an unstable point of equilibrium, for which the exact period diverges. 



VI. CONCLUSIONS 



We have presented a method for the solution of nonlinear problems which are conservative 
and periodic. It is based on the application of the Linear Delta Expansion to the Lindstedt- 
Poincare method. We applied it to two problems: the Duffing Equation and the nonlinear 
pendulum. In the case of the Duffing equation we find that the new model converges faster 



and with greater accuracy than the simple LP method. Also, by comparing it with methods 
based on the perturbative 5 expansion, we show that our solution not only converges faster 
and more accurately, but it also works for a much wider range of parameters, including the 
case in which the nonlinear coupling // is negative. In a similar fashion, we show that the 
method works remarkably well for the solution of the nonlinear pendulum, for which the 
method is implented without performing any Taylor expansion of the potential. 

We are currently workingon the extension of the present method to quantum systems 
and multiple scale analysis 
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APPENDIX A: COEFFICIENTS 

In this Appendix we present the computation of the coefficients of s^+i hi Eq. Let 
us rewrite Eq. (|62|) in the following form: 




(Al) 



We now proceed to compute each of these terms: 



2n+l 




(A2) 



+ 




'2n+l 



+ a 3 A 5 n0 



A 4 A 2 



(A3) 




j(B) 
£ 2n+l 



2 oo oo 

X! X! ^ ^2j+l [<Wi-j,0 + $n+l-l+j,0 + <Wj-n,0 + $n+l+j+l,0 
1=0 j=0 



,-(2) 

~ A 4 A 2 



The four different integrals become: 



2 oo 



T (B1) U V- ~ ,(2) 

'2n+l - ?T C 2(j-n) «2j+l 



2 - 



b) 



,(B2) 
'2n+l 



2 oo 



9 2^ C 2' "2(i-n-l)+l 
Z l=n+l 



~ c) 



r(B3) 
( 2n+l 



u ~ ~ A 2 ) 

2 2^ C 2i »2(n-i)+l 



2 n 



«=0 



- d) 



J 2n+1 — U 



Therefore we finally have that 

,2 



t(B) 
L 2n+1 



c 2(j-n) 4j+l + H C 2« 4(Z-n-l)+l + H C 2« d 2(n-l)+l 
j=n l=n+l 1=0 



7(2) 



/■(C) 
f 2n+l 



/■(CO 

'2n+l 



1 p2tt f oo \ 

= £ / q e* ( 2 "+^ A 2 ( £ cos(2j + l)rj 



2 J(2) 



A d. 



2n+l 



• I {D) 

• J 2n+1 



(D) 1 & 



hn+l 

7T 



r2n ;2 / oo \ 

^ rfr e* {2n+1)T y sin A cost f £ dg+i cos ( 2 -? + X ) r j 

/ rfr e * (2n+Dr £ C2/+i cos(2/ + l)r (AU) 



=0 

00 00 

x E 4!i+i cos(2m+l)r]T d$ +1 cos(2j + l)r 



m=0 j=0 

We need to calculate the following integral: 
1 f 2lT 



J = - / rfre l(2 " +1)r cos[(2j + l)r] cos [(21 + l)r] cos[(2m+l)r] . (A12) 

7T JO 

Using the relation 



C = cos[(2j + l)r] cos[(2/ + l)r] cos [(2m + l)r] 

= i [cos(2(Z + m + j) + 3)r + cos(2(/ + m - j) + l)r 

+ cos(2(/ -m + j) + l)r + cos(2(/ - m - j) - l)r] , (A13) 



one obtains 

1 /-27T 1 , 

2" = - J o dT e %{ - 11 > T C = - |52(rt+Z+m+j)+4,0 + ^2(n-J-m-j)-2,0 
+ ^2(n+«+m-i)+2,0 + ^2(n-/-m+j),0 + ^2(n+/-m+j)+2,0 + ^2(n-/+m-j),0 
+ ^2(n+Z-m-j),0 + &2(n-l+m+j)+2fi} , (A14) 

and finally 

^J2 00 00 00 

4n+l — "7T XI X! X! C2 «+l ^2m+l {^2(n+«+m+j)+4,0 + ^2(n-/-m-j)-2,0 

8 i=0m=0j=0 

+ ^2(n+i+m-j)+2,0 + ^2( rt— i— m+i),0 + ^2(n+«-m+j)+2,0 + $2(n-l+m-j),0 

— ^4 • (A15) 

We are then left with 8 integrals that can be evaluated in the following way (we call 
them Ijj): 

- i) 



^2 00 00 00 

Id — "7T X] X! X! C2 '+! ^2m+l ^2j+l ^2(n+i+m+i)+4,0 = (A16) 
Z=0 m=0 j=0 
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iii) 



iv) 



v) 



- vi) 



r (6) 
l D 



,2 oo oo oo 



^ l_^J LaJ 

r (2) _ ^ Y^ Y^ Y^ j(1) j(1) r 

J D — "7T 2^ 2^ 2^ C2 '+l "2m+l «2j+l °2(n-l-m-j)-2,0 

° 1=0 m=0 j=0 

^2 n-j-ln—1 

= "IT X! X! C 2(n-m-j-l)+l ^2m+l ^2j+l ( A l 7 ) 

m=0 j=0 



. ,2 oo oo oo 

r (3) u j(1) j(1) r 

J D — 2^ 2^ 2^ C 2'+! a 2m+l "2j+l °2(n+«+m-j)+2,0 



Q 

° «=0 m=0 j=0 



^2 oo oo 

"g" X! c 2(j-n-m-l)+l d 2rn+1 d 2j+ i (A18) 

m=0 j=n+m+l 



, ,2 oo oo oo 

(4) _ W y^ ^( : ) A 

J D — ^~ 2^ 2^ 2^ C2/ +! a 2m+l «2j+l °2(n-«-m+j),0 



° /=0m=0jr'=0 
^2 oo n+i 

"7T C 2(n-m+j)+l d 2m+1 d 2j+l (A19) 

j=0m=0 



2 oo oo oo 

J '2m+1 a 2j+l °2(n+l-m+j)+2,0 



— "7T Y Y C 2«+l 4m+l d2j+l $ 



o 

° Z=0 m=0 j=0 

, ,2 oo oc 

W — — ,(1) ,(1) 



— J2 c 2(m ^j_ 1)+l 4/+1 4m+i (A20) 

m=n+j+l j=0 



^jl oo oo oo 

"7T X! X! X! C 2'+l <^2m+l ^2j+l ^2(n-«+m-j),0 

Z=0 m=0 jr=0 

^2 oo m+n 

~g~ Y Y C 2(m+n-j)+l d 2j+1 d 2m+1 (A21) 



m=0 j"=0 



vii) 



, ,2 oo oo oo 

r (7) _ u Y^ Y^ Y^ JX) x 

J D — "7T 2^ Z^ Z^ C 2/+l «2m+l «2j+l °2(n+i-m-j),0 



Q 

° l=0m=0i=0 
. ,2 oo 

W — — 7 (1) j(1) 



X! X! C 2(m+j-n)+l ^2j+l ^2m+l (A22) 

m=max(0,n—j) j=0 



— viii) 



2 oo oo oo 

r (8) ^ Y^ Y^ Y^ ^C 1 ) r 

1 D — ~~cT Z^ Z^ C 2'+! "2m+l «2j+l °2(n-«+m+j)+2,0 

Z=0 m=0 j=0 

^2 oo oo 

— "7T X! X C 2(n+m+j+l)+l ^2j+l ^2m+l (A23) 
8 n o— n 



m=0 j=0 

The final expression is: 

2 [n-j-ln-1 oo oo 

r (D) ^ J Y" Y^ ^C 1 ) Y^ Y" ,/W ^C 1 ) 

J 2n+1 — "3" S Z^ C 2(n-m-j-l)+l «2m+l Q 2j+1 + Z^ Z^ c 2(j-n-m-l)+l «2m+l «2j+l 

8 [ m=0 j=0 m=0 j=n+m+l 

oo n+i oo oo 



+ X] X] C 2(n-m+i)+l ^2m+l ^2j+l + X X C 2(m-n-j-l)+l ^2j+l <^2m+l 
j=0 m=0 m=n+j'+l j=0 

oo rre+n oo oo 

C2(m+n-j)+l ^2j+l ^2m+l C2(m+j-n)+l <^2j+l <^2m+l 

m=0 j=0 m=max(0,n-j) j=0 

oo oo 



+ X X C 2(«+m+j+l)+l 4j+l 4m+l f (A24) 
m=0 i=0 I 



[1] A. Lindstedt, Mem. de l'Ac. Imper. de St. Petersburg 31, 1883 

[2] C.M.Bender, K.A.Milton, S.S.Pinsky and L.M.Simmons, J.Math.Phys.30 (7), 1989 

[3] CM. Bender and L.M.A. Bettencourt, Phys. Rev. D 54, 7710 (1996) |arXiv:hep-tF/ 9607074 . 

[4] A. Okopiriska, Phys. Rev. D 35, 1835 (1987); A. Duncan and M. Moshe, Phys. Lett. B 215, 
352 (1988) 

[5] M. P. Blencowe and A. P. Korte, Phys. Rev. B 56, 9422 (1997) |arXiv:cond-mat79 706260 . 

[6] J. L. Kneur, M. B. Pinto and R. O. Ramos, |arXiv:cond-mat7 0207295 , 

[7] J. L. Kneur, M. B. Pinto and R. O. Ramos, Phys. Rev. Lett. 89, 210403 (2002) 
|arXiv:cond-mat/0207089| . 

[8] G. Krein, D. P. Menezes and M. B. Pinto, Phys. Lett. B 370, 5 (1996) |arXiv:nucl-th/ 9510059 . 

[9] M. B. Pinto and R. O. Ramos, Phys. Rev. D 60, 105005 (1999) |arXiv:hep-ph/9903353| . 

[10] P. Amore and A. Aranda, submitted to Phys. Lett. A (2003) [arXiv:mat£ph /0303042 

[11] P. Amore and A. Aranda, work in progress 



